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摘要 : 脉冲 星 信号 在 星际 空间 传播 的 过 程 中 , 由 于 星际 介质 的 存在 造成 观测 到 的 脉冲 星 
信和 号 发 生 色散 , 因此 需要 对 接收 的 脉冲 星 信号 进行 消 和 色散, 以 获得 原始 的 脉冲 星 信号 。 目 前 ， 
消 色 散 方法 主要 分 为 两 种 ， 相 干 消 色 散 和 非 相 干 消 色 散 。 相 对 来 说 相干 消 色 散 效果 彻底 ， 算 
法 较 简 单 ,， 而 且 能 保留 原始 数据 的 时 间 分 准 滨 ,不 过 计算 量 较 大 , 但 是 现在 快速 进步 的 计算 
i 径 使 计算 量 的 问题 得 到 很 好 的 解决 。 为 了 精确 了 解 两 种 消 色散 方法 的 区 别 , 利用 相 

关系 数 的 方法 定量 地 比较 了 相干 消 色 散 、 非 相干 消 色 散 两 种 方法 的 效果 : 在 一 定 的 频率 值 之 
下 ， 前 者 得 到 消 色 散 效 果 优 于 后 者 。 同 时 确定 两 种 消 色 散 方法 在 效果 相同 时 的 观测 频率 。 
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脉冲 星 是 一 种 具有 超 高 温 、 超 高 压 、 超 高 密度 、 超 强 磁场 、 超 强 电场 和 超 强 引力 场 等 极 
端 物理 条 件 的 天 体 ， 一 般 认 为 是 超新星 爆发 的 产物 ， 典 型 半径 约 为 10 km， 而 质量 却 与 太 P 
相当 , 核心 密度 高 达 1014 g .cm-30。 目 前 , 大 多 数 人 认为 , 脉冲 星 是 高 速 旋转 的 磁 中 子 星 ， 
其 自转 轴 与 磁 轴 之 间 有 一 个 夹 角 , 两 个 磁极 各 有 一 个 辐射 束 。 脉冲 星 的 辐射 几乎 占据 了 整个 
电磁 波谱 ,包括 射电 、 红 外 、 可 见 光 、 紫 外 、X 射线 和 7Y 射线 等 电磁 波 频 段 。 其 辐射 具有 较 
高 的 周期 性 ， 辐 射 周期 和 自转 周期 一 致 ， 且 周期 变化 率 较 小 。 
脉冲 星 极端 的 物理 特性 和 辐射 的 高 周期 性 具有 较 大 研究 意义 和 应 用 前 景 所 。 从 科学 研究 
来 讲 , 脉冲 星 可 以 用 于 验证 和 发 展 极 端 条 件 下 的 物理 学 和 天 文学 理论 , 毫秒 脉冲 双星 的 研究 
己 经 检验 了 广义 相对 论 和 引力 波 辐射 1H， 从 工程 应 用 来 讲 ， 毫 秒 脉冲 星 可 以 用 于 精确 计时 
和 星际 导航 WI。 由 于 毫秒 脉冲 星 的 周期 变化 率 小 于 10-19 s. s-!， 周 期 噪声 极 小 ， 周 期 稳定 
度 比 较 高 ， 且 随时 间 的 增加 会 持续 改善 ,甚至 能 超过 氧 原子 钟 ， 多 个 不 同 天 区 的 毫秒 脉冲 星 
计时 观测 可 以 组 合成 一 个 平均 的 脉冲 星 钟 ， 这 可 以 用 于 现 有 原子 钟 长 期 稳定 性 的 修正 。 
为 了 进行 上 述 研究 ， 就 要 获得 精确 的 脉冲 星 脉 冲 到 达 时 间 中 (Time Of Arrival, TOA), 
由 于 星际 介质 的 存在 ,观测 设备 接收 到 的 脉冲 星 信 号 发 生 色 散 。 为 了 简化 理论 , 星际 介质 被 
看 成 是 低温 、 低 密度 等 离子 体 ， 它 对 在 其 中 传播 的 频率 为 1 的 电磁 波 的 折射 率 为 
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对 于 射电 波段 观测 , 接收 到 的 脉冲 信号 是 一 维 时 域 上 的 脉冲 轮廓 , 不 同 频率 的 脉冲 之 间 
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其 中 , DM (Dispersion Measure) 为 传播 路 径 上 的 星际 介质 的 色散 量 , 色散 量 DM = i ne dl = 


元 dd， 单位 为 cm 3 .pc 
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元 为 传播 路 径 上 的 平均 电子 数 密度 ; D 为 星际 介质 的 色散 常量 ，D = 一 一 = 


4148.808MHz2pc-tcm3s; 

fiow 和 hign 分 别 为 接收 频带 的 起 始 频 率 和 截止 频率 ， 单 位 为 MHz。 
(2) 式 可 知 ， 一 定 带 宽 脉 冲 星 信 号 的 色散 延迟 会 造成 时 域 脉冲 的 相位 移动 、 幅 度 降低 和 脉 
冲 轮廓 展 宽 等 色散 现象 , 这 限制 了 到 达 时 间 分 辨 率 。 因 此 需要 对 接收 到 的 脉冲 星 信 号 进行 消 
色散 并 获得 精确 的 到 达 时 间 ， 通 过 对 脉冲 到 达 时 间 的 拟 合 可 以 得 到 脉冲 星 的 相关 物理 参数 ， 
进而 研究 其 自身 的 辐射 机 制 和 其 他 相关 特性 。 


jami 


N 通道 滤波 器 组 


图 1 消 色 散 算法 流程 图 ， 上 半 部 分 是 相干 消 色 散 的 流程 图 ， 下 半 部 分 是 非 相干 消 色散 的 流程 图 


Fig.1 Flow chart of coherent de-dispersion algorithm; the above belongs to coherent de-dispersion, the below belongs to incoherent 


de-dispersion 


消 色 散 的 宗 则 是 消除 信号 频带 内 的 色散 效应 ,目前 主要 有 非 相 干 消 色 散 和 相干 消 色散 叫 。 
非 相 干 消 色 散 是 基于 时 间 序 列 , 用 滤波 器 组 分 通道 ,， 按 相 对 于 参考 通道 的 延迟 时 间 单 个 通道 
进行 平移 从 而 消除 色散 ， 如 图 1， 相 干 消 色 散 是 在 频 域 上 利用 星际 介质 的 传输 函数 (chirp R 
数 ) 实 现 通 带 内 的 消 色散 ， 如 图 1。 
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图 2 脉冲 星 PSR J2322+2057 的 脉冲 轮廓 。 实 线 是 相干 消 色 散 后 的 轮廓 ， 虚 线 是 非 相 干 消 色 散 后 的 轮廓 外 
Fig.2 Pulse profile of PSR J2322+2057. Solid line: pulse profile observed by coherent de-dispersion. Dashed line: The same profile, 


filtered by incoherent de-dispersion 


目前 已 经 有 学 者 在 比较 两 种 消 色 散 方法 方面 进行 了 一 些 研究 。 文 [6] 在 430MHz 处 对 PSR 
J2322+2057 进行 了 观测 、 获 得 了 相干 和 非 相 干 消 色 散 后 的 平均 脉冲 轮廓， 并 计算 出 脉冲 到 
达 时 间 的 精度 ， 后 者 与 前 者 到 达 时 间 的 均 方 根 误差 比值 为 0.3 ， 如 图 2。 文 [7] 用 DSP 算法 
在 327MHz 处 比较 了 这 两 种 方法 得 到 的 一 些 脉 冲 星 信 号 的 积分 轮廓 ， 相 干 消 色 散 后 数据 的 
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时 域 序列 脉冲 幅度 和 时 间 分 辨 率 是 非 相 干 消 色 散 后 数据 的 10~15 倍 中 。 国 内 ， 中 国 科学 院 
新 疆 天 文 台 的 脉冲 星 观测 团队 利用 基于 MK5A 的 VLBI 记录 系统 建立 一 套 带宽 为 64MHz 的 
相干 消 色 散 系统 , 得 到 了 一 些 脉冲 星 的 平均 脉冲 轮廓 ， 与 以 前 使 用 的 2x128x2.5MHz 多 通道 
非 相 干 消 色 散 系统 的 结果 进行 了 比较 外， 结论 是 相干 消 色 散 优 于 非 相 干 消 色 散 ， 但 没有 进行 
系统 和 定量 的 分 析 研 究 。 本 文 旨 在 用 相关 系数 的 方法 定量 地 分 析 两 种 方法 消 色散 效果 的 差别 。 


1 消 色 散 原 理 

两 种 消 色 散 方法 基于 不 同 的 消 色散 原理 。 非 相干 消 色 散 是 基于 时 域 上 各 个 频段 的 信号 相 
位 移动 进行 的 不 连续 消 色散 ， 相 干 消 色 散 是 基于 频 域 乘积 的 连续 消 色 散 ， 具 体 描述 如 下 。 
1.1 非 相干 消 色 散 

非 相 干 消 色 散 的 过 程 是 首先 把 观测 带宽 为 BW 的 时 域 信号 由 滤波 器 组 分 成 若干 狭窄 通 
道 , 然后 在 时 域 上 把 每 个 通道 信号 按 适当 的 时 延 量 进行 平移 , 最 后 将 脉冲 对 齐 后 的 所 有 通道 
言 号 幅度 累加 得 到 消 色 散 后 的 时 域 信号 序列 ， 如 图 3。 非 相干 消 色 散 的 优点 是 数据 量 小 ， 计 
算 量 少 。 这 种 方法 的 缺点 : (1) 单个 通道 内 的 色散 效应 没有 消除 ， 增 大 通道 的 数量 可 以 减少 
剩余 的 色散 效应 ， 但 是 由 不 确定 性 原理 tes。 x 1/(2 x bw) 可 知 ， 增 大 通道 数量 ， 单 个 通道 信 
号 的 时 间 分 辨 率 会 增 大 ， 所 以 要 根据 需求 ,适当 调整 通道 数 的 数值 ，(2〉 由 于 非 相 干 消 色散 
过 程 是 单 通道 时 域 序 列 的 平移 ， 该 方法 会 损失 信和 号 的 相位 信息 。 
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图 3 消 色 散 原理 图 四。 上 半 部 分 显示 的 是 非 相干 消 色 散 的 原理 ， 下 半 部 分 显示 的 是 相干 消 色 散 的 原理 


Fig.3 Schematic diagram of incoherent de-dispersion. The upper part shows the principle of incoherent de-dispersion, the utter part shows 


the Principle of coherent de-dispersion. 


1.2 相干 消 色 散 
星际 介质 的 色散 效应 相当 于 移 相 器 的 作用 , 因此 接收 到 的 脉冲 星 信 号 相当 于 原始 信号 加 
上 移 相 器 得 到 的 结果 ， 移 相 器 可 以 用 传递 函数 只 (六 来 表征 ， 有 具体 表达 中 为 
i2nD:DM-f 47 
H(f, + fo) =e Attori (3) 
KE, fA; 万 为 中 频频 率 。 
由 传递 函数 瓦 (月 可 知 ， 原 始 信 号 可 以 通过 乘 以 其 复 共 箔 号 (一 用 进行 完全 消 色 散 ， 消 除 
整个 观测 带宽 内 的 色散 相应 ， 如 图 3， 最 大 的 时 间 分 辨 率 为 1/(2. BW)。 这 种 消 色 散 方 法 可 
用 于 高 精度 毫秒 脉冲 星 计 时 和 脉冲 星 射 电 辐 射 过 程 的 研究 。 
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理论 上 讲 , 两 种 消 色 散 结 果 的 差异 主要 由 各 个 通道 的 剩余 色散 量 造 成 的 累加 效应 , 因此 
分 析 两 种 消 色 散 效果 可 以 通过 分 析 各 个 通道 的 剩余 色散 量 实现 。 
对 于 上 边 带 信号 的 非 相干 消 色 散 ， 由 (2) 式 可 知 ， 整 个 带宽 BW 内 的 最 大 时 延 : 


2foọ+BW):BW 
tpeiay = D* DM- Cer te (4) 
单个 通 EaR 3 
2:bw 
t = D:DM- ， 5 
“i Motnow e ow _ 
其 中 ，bw 为 单个 通道 带宽 ，bw = So N_binf 为 滤波 器 组 的 通道 数 ，n 的 取 值 范围 为 [1， 
N_binf /2]. 对 于 fo > bw, 
2:bw 
faery ~ DDM Trapwp 二 


2 脉冲 星 信号 模拟 、 消 色散 算法 及 比较 方法 实现 


本 研究 用 MATLAB 程序 模拟 了 被 射电 望远镜 接收 到 的 脉冲 星 信号 ， 实 现 了 两 种 消 色散 
过 程 和 两 种 消 色 散 方法 效果 的 定量 比较 过 程 ， 有 具体 过 程 如 下 。 
2.1 原始 信号 模拟 
脉冲 星 辐 射 的 原始 信号 很 弱 ， 为 了 提高 信 噪 比 ， 需 要 长 时 间 的 积分 。 为 了 简化 计算 ,此 
处 模拟 一 个 周期 的 基带 信号 。 为 了 便于 比较 两 种 方法 的 消 色 散 效 果 ， 周 期 采用 短 周期 
0.00390625 秒 ， 同 时 为 了 提高 程序 的 运算 效率 ,采样 率 (2Gbps)、 采 样 时 长 和 信和 号 带宽 BW 
等 都 取 2 的 容 指 数 。 加 色散 前 的 基带 信号 和 加 色散 后 基带 信号 如 图 4。 
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(a) (b) 
图 4 (a) 加 色散 前 的 基带 模拟 信号 ，(b) 加 色散 后 的 基带 模拟 信号 ; 信号 频谱 带宽 BW 为 1GHz， 起 始 频率 太 为 1GHz， 色 散 量 
DM 为 40 cm-3 .pc， 图 5 也 取 同 样 的 取 值 


Fig.4 The left image is the baseband analog signal before adding dispersion, the right image is the baseband signal after the adding 


dispersion; Signal spectrum bandwidth BW : 1 GHz, the starting frequency f, : 1 GHz, the Dispersion Measure DM : 40 cm-3 - pc. In the 


figure. 5 those variables also take same value 


(1) 模拟 带宽 为 BW, ALIEN fo, 谱 指数 为 -2.0 的 基带 信和 号 口 。 
(2) 把 频 域 的 模拟 基带 信号 通过 逆 傅 里 叶 变 换 到 时 域 ， 进 行 高 斯 脉冲 调 秆 
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(3) FE ved till Jes fies SFR a ee I Ae St, EW aT eR A H( 有 后 进行 加 色散 。 
(4) 把 加 色散 的 频 域 信号 再 变换 到 时 域 ， 加 高 斯 白 噪声 成 模拟 的 脉冲 星 时 域 信号 。 
2. 2 消 色散 算法 实现 
根据 模拟 长 度 为 N_data 的 脉冲 星 信 号 数据 , 利用 MATLAB 程序 实现 信号 的 消 色 散 过 程 ， 

并 实现 两 种 消 色 散 后 的 序列 与 加 色散 前 序列 的 相关 过 程 。 
2. 2. 1 非 相 干 消 色散 步骤 

(数据 分 箱 : 把 模拟 的 基带 数据 按 顺序 平均 分 箱 。 

(2) 多 相位 滤波 010: 对 每 个 箱 中 的 数据 进行 多 相位 滤波 ， 数 据 序列 长 度 减 小 原始 


的 2; 这 可 以 减弱 传 里 叶 变 换 的 频谱 泄露 。 多 相位 滤波 过 程 为 : 利用 相同 长 度 的 sinc 
函数 序列 依次 乘 以 每 个 箱 中 的 数据 , 把 得 到 的 序列 按 顺 序 平均 分 成 4 段 ， 并 逐次 对 
应 相 加 得 到 原始 箱 内 序列 :长 度数 据 。 

(3) 傅 里 叶 变换 ， 对 每 个 箱 中 的 数据 进行 傅 里 叶 变 换 。 

(4) 分 通道 并 平移 数据 : 取 参 考 频率 为 截止 频率 刻 ion， 任意 通道 中 心 频率 /的 通道 按 
延迟 时 间 平 移 。 由 (2) 式 可 知 ， 相 对 延迟 时 间 为 


1 1 

P i Á 
(5) 逆 传 里 叶 变换 并 求 和 : EAA PEET ARR, RA A AA 
据 求 和 得 到 分 箱 后 的 数据 序列 ， 非 相干 消 色 散 后 的 信号 见 图 5。 
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图 5 (a) 非 相干 消 色 散 后 信号 的 二 维 频谱 ; (b) 相 干 消 色散 后 信号 的 二 维 频谱 


Fig.5 The two-dimensional spectrum of the signal after incoherent de-dispersion, the two-dimensional spectrum of the signal after 


coherent de-dispersion 


2. 2. 2 相干 消 色 散步 又 
基带 信号 只 有 一 个 周期 , 为 了 简化 运算 ,将 每 次 处 理 的 数据 长 度 也 定 为 一 个 周期 。 具体 
步骤 如 下 
(直接 进行 傅 里 叶 变 换 : 把 模拟 的 基带 数据 进行 傅 里 叶 变 换 ， 得 到 频 域 数据 。 
(2) 频 域 消 色散 : 把 频 域 数据 乘 以 消 色 散 函 数 H( 一 站 ) 得 到 消 色散 后 的 频 域 数 据 。 
(3) 进 行 逆 傅 里 叶 变 换 : 把 消 色散 后 的 频 域 数 据 进行 逆 傅 里 叶 变 换 得 到 原始 序列 长 度 的 消 色 
散 后 的 时 域 数据 。 
(Apr eM a PI: 把 消 色散 后 的 数据 依次 折 麦 成 与 非 相干 消 色 散 分 箱 后 序列 长 度 相同 的 


箱 ， 然 后 把 每 个 箱 中 的 数据 累加 得 到 序列 ， 相 干 消 色散 后 的 信号 见 图 5。 
2. 3 消 色散 效果 的 定量 比较 

为 了 定量 分 析 两 种 消 色 散 , 让 两 种 消 色 散 后 的 结果 与 加 色散 前 的 数据 进行 互相 关 , RE 
们 的 互相 关系 数 : 


wot 1Xi'Vi) 


For ELA 


其 中 ，xi 属 于 长 度 为 INS X, yJRFRENN 的 序列 Y。 目 前 普遍 认为 ， 相 关系 数 的 绝 
对 值 在 0.3 以 下 代表 无 直线 相关 ，0.3 以 上 代表 直线 相关 ，0.3-0.5 代表 低 度 相 关 ，0.5-0.8 代 
表 显 著 相 关中 等 程度 相关 )，0.8 以 上 代表 高 度 相 关 ; 当 且 仅 当 其 绝对 值 为 1 时 ,代表 线性 
相关 。 对 于 研究 的 信号 ， 表 明 两 个 信号 完全 相同 ， 如 相同 的 脉 宽 和 信 噪 比 。 有 具体 步骤 为 ; 
(1) 把 加 色散 前 的 模拟 数据 按 上 面 的 方法 折 又 累加 成 长 度 与 非 相 干 消 色 散 分 箱 后 序列 长 度 相 
同 的 数据 序列 。 
(2) 把 两 种 消 色 散 后 的 数据 与 上 一 步 的 分 箱 数 据 进 行 互 相关 得 到 互相 关系 数 。 
3 两 种 消 色 散 的 定量 比较 
首先 ， 根 据 傅 里 时 变换 的 理论 ,给 出 两 种 消 色 散 过 程 的 主要 计算 量 ， 并 记录 两 种 消 色散 
过 程 的 耗 时 ; 其次， 定量 比较 两 种 消 色散 的 效果 ; 最 后 ， 两 种 消 色 散 效 果 相 当时 ， 估 计 信 和 号 
频谱 的 起 始 频 率 。 
在 信号 模拟 和 消 色 散 过 程 中 ,色散 量 的 取 值 为 1~40 cm- pce, HEA cm: pe. fa 
号 频谱 带宽 BW 取 值 为 32 MHz, 64 MHz、128 MHz、256 MHz、512 MHz 和 1 GHz， 单 通 
道 带宽 为 1 MHz。 信 号 频谱 起 始 频率 的 取 值 范围 为 236 MHz ~ 8 GHz, 对 于 比较 运算 耗 时 和 
消 色 散 效果 , 其 步 长 为 64 MHz， 对 于 确定 两 种 消 色散 效果 相同 时 的 频谱 起 始 频率 , 其 步 长 为 
16 MHz 。 
3.1 两 种 消 色 散 结 果 得 到 的 互相 关系 数 的 比较 
通过 计算 得 到 两 种 消 色 散 后 的 数据 与 加 色散 前 数据 的 相关 系数 如 图 6、 图 7。 首先 分 析 
两 种 消 色 散 结果 得 到 互相 关系 数 随 各 变量 的 变化 特点 。 
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q6 相干 消 色散 后 的 数据 与 加 色散 前 的 数据 的 相关 系数 ， 为 了 便于 显示 ， 此 处 相关 系数 乘 了 1000 倍 


Fig.6 the coherence coefficient of the data after the coherent de-dispersion with the data before the adding dispersion. To facilitate the 


display, the coherence coefficient is multiplied by 1000 


由 图 6 可 知 ， 相 干 消 色散 得 到 的 相关 系数 接近 于 1， 这 说 明 ， 在 不 考虑 离散 化 运算 过 程 
中 引入 的 误差 时 ， 相干 消 色散 可 以 完全 恢复 原始 加 色散 之 前 的 信号 。 结 合 两 图 可 以 发 现 ， 当 
BW 减 小 时 ， 相 关系 数 的 波动 变 大 ， 这 与 (9) 式 有 关 ; 其 次 ， 数 据 处 理 过 程 存在 频谱 泄露 ， 这 


是 由 于 数据 在 周期 性 分 箱 过 程 中 产生 一 些 截 断 频 点 , 有 限 长 离散 数据 的 传 里 叶 变 换 有 频谱 污 
露 ， 这 些 效应 的 影响 随 着 BW 的 减 小 而 增 大 。 

两 种 消 色散 得 到 的 相关 系数 见 图 7。 对 于 相同 的 色散 量 PDM， 频 谱 起 始 频率 万 的 取 值 在 
一 定 范围 内 ， 相 关系 数 之 差 随 着 变量 fp 指数 减 小 ， 但 到 某 值 户 之 上 ， 该 相关 系数 之 差 接 近 
于 1 且 波 动 逐渐 变 小 ， 这 说 明 在 频谱 起 始 频率 FKP 有 时 ， 非 相干 消 色 散 的 效果 比 相干 消 
色散 的 差 ， 且 两 种 方法 的 消 色 散 效 果 的 差距 随 着 的 增加 而 缩小 ， 在 有 高 于 有 时， 两 种 方 
法 的 消 色 散 效果 相同 。 结 合 两 图 可 以 看 出 ， 相 同 万 、 随 DM 变化 对 应 的 相关 系数 的 变化 曲线 
随 着 BW 的 减 小 起 伏 变 大 ， 这 是 数据 处 理 中 的 频谱 泄露 的 误差 造成 的 。 相 同 有 、 相 同 DM 对 
应 的 相关 系数 之 差 随 BW 的 增 大 而 减 小 , H. fi 随 BW 的 减 小 而 增 大 。 结合 (6) 式 和 图 7 可 知 ， 
在 不 考虑 误差 范围 时 ,在 一 定 观 测 频率 内 ， 相 干 消 色散 效果 优 于 非 相干 消 色散 效果 ， 两 种 方 
法 的 消 色 散 效果 的 差距 随 观 测 带宽 、 观 测 频率 的 增加 而 按 指数 关系 减 小 , 而 随 着 色散 量 的 增 
大 而 按时 指数 关系 增 大 ,同时 能 够 确定 ,可 以 找到 两 种 方法 消 色 散 效果 相同 时 的 观测 起 始 频 
率 。 
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图 7 两 种 消 色 散 方法 得 到 的 相关 系数 之 差 


Fig.7 The coherence coefficient of the data after the incoherent de-dispersion with the data before the adding dispersion 


3.2 两 种 消 色 散 效 果 相 同 对 应 的 通 带 起 始 频率 
两 种 消 色 散 效 果 相同 时 , 这 对 应 于 非 相 干 消 色 散 单 通道 内 的 时 延 taerty 取 某 一 接近 零 常数 。 
非 相 干 消 色 散 过 程 ， 若 通道 带宽 bw 为 常数 ， 由 (6) 式 可 得 


(9) 


oy gee 
dely ` [f +nbw]? 


jami 


HER, taewy 为 常数 ， 有 


fa +nbw x YDM . (10) 

经 过 计算 , 在 两 种 消 色 散 后 的 数据 与 加 色散 前 的 数据 的 相关 系数 小 于 茶 一 定 值 条 件 下 ,得 到 
频谱 起 始 频率 万 如 表 1， 考 虑 到 如 上 所 述 数据 处 理 中 引入 的 误差 ， 这 个 值 取 0.01. 

4 结论 

运用 MATLAB 程序 实现 了 两 种 消 色 散 过 程 的 定量 比较 ， 并 确定 了 两 种 消 色散 方法 相同 

时 的 频谱 起 始 频 率 , 在 实现 比较 过 程 , 色散 量 的 取 值 为 1~40 cm -3 pe, 步 长 为 1 cm. pc; 

信号 频谱 带宽 BW 取 值 为 32 MHz、64 MHz、128 MHz、256 MHz、512 MHz 和 1 GHz, % 

通道 带宽 为 1 MHz; 信号 频谱 起 始 频率 的 取 值 范围 为 256 MHz ~ 8 GHz， 对 于 比较 运算 耗 时 

和 消 色 散 效 果 , 其 步 长 为 64 MHz， 对 于 确定 两 种 消 色 散 效 果 相同 时 的 频谱 起 始 频 率 ,其 步 长 

为 16 MHz。 定 量 的 比较 结果 说 明 在 一 定 的 观测 频率 之 内 相干 消 色 散 的 效果 优 于 非 相 干 消 色 


散 ， 实 际 观测 中 可 以 根据 观测 要 求 ， 参 考 这 些 结果 ， 选 择 消 色散 方法 ; 最 后 研究 了 在 不 同色 
散 量 和 不 同 频谱 带宽 下 , 两 种 消 色 散 技术 消 色 散 效果 相同 时 的 起 始 频率 , 非 相 干 消 色散 脉冲 


星 观测 系统 在 实际 的 观测 中 可 以 参考 这 些 频率 制定 观测 频段 。 非 相干 消 色散 过 程 引 用 了 多 相 


位 滤波 器 ， 一 定 程 度 上 缓解 了 传 里 叶 变换 的 频谱 泄露 问题 ; 但 是 频谱 泄露 问题 依然 存在 ， 尤 


N 


在 低频 端 和 高 频 端 。 下 一 步 将 扩展 色散 量 的 变化 范围 和 变化 的 单个 通道 带宽 , 同时 继续 改 


善 频谱 泄漏 问题 。 


表 1 不 同 BN, OM 对 应 消 色 散 效 果 相 同时 的 A, AB firs MHz; 这 只 是 部 分 数据 


Table 1 the f, to different BW and DM, which corresponds to the de-dispersion effect, the 


unit is MHz; This is just part of the data 


BW / MHz 32 64 128 256 512 1024 
DM=5 cm™?:pc 608 592 560 512 448 400 
DM =10 cm-3.pc 800 736 704 656 592 528 
DM =15 cm-3.pc 912 848 816 768 688 608 
DM =20 cm-3.pc 976 928 896 848 768 688 
DM =25 cm-3.pc 1056 1024 992 912 832 752 
DM =30 cm-3.pc 1152 1072 1056 976 896 800 
DM =35 cm-3.pc 1232 1152 1104 1040 960 848 
DM = 40 cm-3.pc 1232 1200 1152 1088 1008 896 
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Comparative Study between the Coherent De-dispersion and 


the Incoherent De-dispersion of Pulsar Signal 


Huang Yuxiang!?, Wang Min!, Hao Longfei!, Li Zhixuan!, Xu Yonghua! 
(1. Yunnan Observatory of Chinese Academy of Sciences, Yunnan Kunming 650011; 2. University 
of Chinese Academy of Sciences, Beijing 100049) 

Abstract: Spreading through in the interstellar space, pulsar signals, due to the existence of 
the interstellar medium, would be observed dispersively, therefore we need to remove dispersion 
of the received pulsar signals for obtaining the original pulsar signals. The current de-dispersion 
methods are mainly divided into two types, coherent de-dispersion and incoherent de-dispersion. 
Relatively speaking, the former is able to eliminate dispersion completely, which has simpler 
algorithm, and preserves the time resolution of the original signal, but has the larger amount of 
calculation. But the rapid progress of computer technology has already solved the problem of 
calculated amount. In order to accurately understand the difference between the two methods of 
de-dispersion, this paper, adopting the method of coherence coefficient, quantitatively compares 
the effects of coherent de-dispersion and incoherent de-dispersion: below a certain frequency, the 
former is better than the latter in obtaining the de-dispersion effects. Meanwhile, the observed 
frequency of the two de-dispersion methods at the same effect is determined. 


Key words: Pulsar; Coherent de-dispersion; Incoherent de-dispersion 


